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FOREWORD 


The  work  described  In  this  report  was  carried  out  by  the  Aerodynamic 
Research  Department  of  the  Cornell  Aeronautical  Laboratory  under  Air  Force 
Contract  Wo,  AF  33(038)<-l?2l*2,  Expenditure  Order  No.  1*60-31-12-15  SB.  lg„  This 
contract  «u  administered  by  the  Aeronautical  Research  Laboratory,  Wright  Air 
Development  Center p under  the  technical  cognisance  of  Hr.  Lee  S„  Wassemain. 

The  experimental  work  carried  out  under  this  contract  »as  supervised  by 
Hr,  John  0,  Wilder  of  the  Cornell  Aeronautical  Laboratory 


ABSTRACT 


Theoretical  and  experimental  studies  of  the  stability  of  flow  in  air 
induction  systems  with  high  pr assure  recoveries  for  boundary-layer  suction 
have  been  carried  out.  The  boundary-layer  flow  entering  such  systems  nay  have 
the  characteristic  of  increasing  its  total  head  with  increasing  flow  rate  for 
part  of  the  operating  range.  This  characteristic  produces  both  static  and 
dynamic  instability.  The  static  instability  is  evidenced  by  the  appearance 
of  unequal  flows  in  ostensibly  identical  branches  of  a system.  The  dynamic 
instability  occurs  in  the  form  of  regular  oscillations  of  the  flow.  These 
oscillations  take  place  at  the  characteristic  frequencies  of  the  duct  and  plenum 
chamber  system.  In  simple  cases,  the  measured  frequency  shows  good  agreement 
with  that  calculated  on  the  basis  of  the  acoustic  theory  of  the  Helmholtz  reson- 
ator. In  the  experimental  investigation,  it  was  shown  that  the  installation  of 
splitter  plates  in  a system  with  a wide  slot  which  exhibited  both  static  and 
dynamic  instability  eliminated  the  latter  but  not  the  former.  Both  types  of 
instability- could  be  eliminated  by  introduction  of  high  slot  and  duct  losses 
into  the  system  or  by  ingesting  a flow  quantity  into  the  system  sufficient  to 
remove  almost  the  entire  boundary  layer  ahead  of  all  slots.  It  should  be 
pointed  out  that  a model  setup  simulating  an  airfoil  boundary  layer  control 
system  would  not  exhibit  either  of  the  instabilitiese  This  was  thought  to  be 
possible  due  to  a plenum  volume  scale  effect  and/or  to  the  relatively  higher 
losses  caused  by  the  model-scale  effect. 
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IHTRGDuCTIOH 


Instability  of  flow  in  Air  induction  systems  which  take  in  boundary- layer 
air  has  often  been  observed,  vainly  in  connection  with  engine  air  intakes.  In 
such  cases,  the  ingestion  of  boundary-layer  air  is  an  incidental  and,  in  soei 
cases,  undesirable  feature  of  the  system.  It  has  often  been  found  possible  to 
cure  the  instability,  and  at  the  a sue  tins  to  improve  the  efficiency  of  the  in- 
take by  diverting  th?  boundary-layer  air  through  a separate  boundary-layer  bleed 
system.  In  the  esse  of  boundary-layer  control  systems,  however,  practically  the 
entire  mass  of  air  being  sucked  through  the  system  is  boundary-layer  air  and  the 
problem  of  stability  of  the  system  needs  more  careful  study.  Although  the  use  of 
boundary-layer  suction  to  improve  the  lift  and  drag  characteristics  of  sirfoils 
is  an  old  concept,  it  has  recently  received  renewed  attention  (refs,  1,2,3  and  k). 
In  particular,  the  possibility  of  stabilisation  of  the  laminar  boundary  layer  by 
suction  appears  attractive  and  is  being  vigorously  pursued.  Although  there  is 
now  a very  considerable  literature  dealing  with  boundary-layer  auction,  the  de- 
tailed consideration  of  the  stability  of  flow  in  the  induction  system  has  been 
almost  entirely  ignored.  Smith  and  Roberts  (ref,  3)  hare,  however,  pointed  out 
that  there  is  a possibility  of  instability  in  boundary-layer  slots  due  to  the 
fact  that,  for  small  flow  rates  into  the  slot,  the  total  head  of  the  air  enter- 
ing the  system  exhibits  a positive  increase  with  increase  of  volume  flaw it 
is  well  known  that  such  a flow  characteristic  is  associated  with  the  phenomena 
of  compressor  surge  (refs,  5 and  6)  and  supersonic  diffuser  "buss",  As  will  be 
seen  later,  this  characteristic  of  the  Blot  intake  is  a necessary  but  not  suf- 
ficient condition  for  instability.  Other  necessary  conditions  for  instability 
arise  from  the  characteristics  of  the  duct  system.  Further,  the  instability 
encountered  may  be  either  static  or  dynamic  in  nature 5 the  type  is  also  deter- 
mined by  the  characteristics  of  the  entire  induction  system, 

Aeith  and  Roberts  (ref0  3)  determined  the  Curve  of  total,  pressure  recovery 
versus  rate  of  flow  for  the  flow  of  boundary  layer  air  into  a slot,  including 
the  estimated  effects  of  slot  losses,  which  are  stabilisiiigo  Their  results 
showed  that  the  slope  of  the  total  pressure  versus  rate-- of- flow  curve  could  be 
positive  for  the  lower  rates  of  flow,  which  could  be  typical  of  those  employed 
for  stabilisation  of  the  laminar  boundary  layer.  Although  the  type  of  instab- 
ility with  which  they  were  concerned  is  not  clear  from  their  paper.  Smith,  in 
a private  communication  to  a member  of  the  C,A0L,  staff  (ref, 11)  seated  that 
the  investigation  was  concerned  with  the  problem  of  long  slots  which  often 
•show  a great  rush  of  air  in  one  side  and  very  slow,  or  even  reverse,  flow  on 
the  other  side".  As  will  be  seen  later,  serious  effects  of  this  type  do  occur, 
but  they  may  also  be  accompanied  by  a resonant  oscillation  of  the  system  as  well. 

The  theoretical  analyses  presented  in  this  report  show  that  the  physical 
factors  underlying  the  static  Instability  dio^ussed  above  may  tuou  give  rise  to 
a dynamic  oscillatory  instability.  It  is  shown  that  the  positive  slope  of  the 
curve  of  total  head  versus  flow  rate  for  the  intake  slot  is  equivalent  to  nega- 
tive damping  imposed  on  the  resonant  system  made  up  of  the  air  in  the  duct  system. 
The  resonant  character  of  this  system  is  vary  much  like  that  of  the  classical 
Helmholts  resonator  in  acoustics  (ref0  8)0  In  fact,  it  was  found  possible  to 
calculate  the  resonant  frequencies  for  the  cases  tested  with  fair  accuracy  based 
on  acoustic  theory. 
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As  far  «3  was  known,  there  ware  no  detailed  experimental  data  on  either  the 
static  or  dynamic  stability  of  boundary-layer  air  induction  systems  available  in 
ths  Aterature.  Since  the  details  of  the  flow  into  a slot  and  particularly  the 
losses  under  dynamic  conditions  cannot  be  predicted  theoretically  with  satisfac- 
tory accuracy^  it  was  difficult  to  ass sea  the  probability  of  encountering  the 
various  types  of  instability  prior  to  conducting  the  present  experimental  program. 
For  this  reason;  much  of  ths  present  experimental  work  was  exploratory  in  nature 
with  the  aim  of  establishing  the  essential  features  of  the  instabilities  encountered 
and  the  mein  conditions  for  their  occurrence.  Also,  the  models  employed  were  Ideal- 
ised representations  of  boundary-layer  control  systems  rather  than  representative 
of  any  practical  case. 

In  the  original  plan  for  the  present  program,  it  was  Intended  that  the  inves- 
tigation be  centered  on  the  question  of  static  instability,  with  only  secondary 
attention  to  the  dynamic  problem,  since  the  existence  and  range  of  occurrence  of 
the  latter  were  uncertain.  Actually,  however,  the  dynamic  instability  presented 
itself  in  the  very  first  tests,  and  it  was  necessary  to  give  primary  attention  to 
it  to  the  detriment  of  the  static  investigations , Pure  static  instability  was 
observed  in  the  present  program  only  with  multiple  air  intakes,  although  a non- 
uniform  distribution  of  the  mean  flow,  of  the  type  described  by  Smith  (ref.  11} , 
was  also  present  during  resonant  oscillation  of  the  system  with  long  narrow  dots. 

In  the  original  program,  it  was  planned  that  most  of  the  work  would  be  done 
with  laminar  boundary  layers  representative  of  boundary-layer  control  for  drag 
r*duetion>.  Tt  turned  out,  however,  that  with  the  small  thickness  of  natural 
laminar  boundary  layer  available  on  an  airfoil  in  these  experiments,  the  losses 
in  the  very  snail  slots  required  may  have  produced  a stable  total  head  versos 
flow  rate  curve,  ioe0,  no  increase  in  total  head  with  increasing  flow  rate.  Thus, 
neither  static  nor  dynamic  instability  could  be  produced  with  a laminar  boundary 
layer  on  the  small  airfoil  model.  It  was  verified  that  a turbulent  boundary  layer 
on  the  same  model  also  gave  no  instability.  The  results  for  instability  reported 
herein  are  therefore  entirely  for  high  Reynolds  number  turbulent  flow  ahead  of 
the  slot.  Such  tests  were  carried  out  by  sucking  the  wall  boundary  layer  from 
the  16*  wall  of  the  C.A.L.  3"  x 16"  subsonic  wind  tunnelo 
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THEORETICAL  ANALYSIS 


Bound ary- Layer  Characteristics 


For  purposes  of  this  analysis,  it  will  ba  adequate  to  ecstnaa  that  the 
boundary-layer  Telocity  profile  is  given  by  a power  law,  naaely 


£-(W 


- (1) 


where  S.  is  the  boundary-vlayer  ^thickness0 « If  the  flow  entering  the  slot 
is  saaller  than  that  corro&ponuing  to  the  entire  botmdary-layer  thickness, 
then  q 9 the  rolime  flow  ,«atesteg  the  slot  trait  tto,  is  given  by 

f 6 r ^ r \ n 

q - ad^  -»  J 'i»  (%  ) dy 

where  8 is  the  thickness  of  flow  id  thin  the  boundary  layer  which  eaters  the 

slot.  Thus, 


**  (hT 


(2) 


The  moan  total  head  of  the  air  entering  the  slot  is  given  by 

,s 


h>  • 4£  ■ u ^ ^ d«# 


where  R is  the  static  pressure  of  the  flew  near  the  slot  or 

8 


H - R 


I fu< 


dy 


low*- 


[ J {B/k)  dy 


■£&  s. 


3n*l 


(3) 


Then 


»n  + i / 


H L * R + 


£f*  ( %.)  / 2 tS^  (%$" 
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or 


H.—  R «■  (%.) 


an 


- (U) 


From  Bq.  (2) 


S/a 


[wj 


( 

r»*T 


Substituting  this  result  into  Bq.  (U),  there  results 


Ho—  R + 


--  (5) 


From  this  equation,  the  Talus  of  *5-tL  which  is  important  in  considerations  of 
stability  may  be  obtained  *r  dq 


or  in  terns  of 


n- 1 


• <a  ® (6) 

• (7) 


If  a laminar  boundary  layer  is  approximated  by  h=»l  , the  value  of  "Sq  is  a 
constant  equal  to  for  all  flow  rates,  as  would  be  expected  for  a 

straight  line  profile.  For  a turbulent  boundary  layer  (assuming  a 1/7  power 
law),  the  result  is 


(8) 


which  slowly  decreases  with  increasing  5 or  q , and  is  smaller  than  the  lam- 
inar value  for  the  higher  flow  rates. 

As  has  been  indicated  in  the  introduction,  instability  is  associated  with 
positive  value  of  . It  would,  however,  be  a mistake  to  conclude  from  J£qs. 

(7)  and  (8)  that  a flow  with  laminar  profile  would  be  more  unstable  than  one  for 
a turbulent  profile,  say  with  $£#»-  , because  the  effects  of  slot  losses  have 

not  yet  been  taken  into  account.  Such  losses,  which  have  to  do  with  the  mixing 
and  diffusion  of  the  boundary  layer  air  on  and  after  entry  into  the  slot,  are 
not  well  understood  and  it  appears  entirely  possible  that  the  losses  for  a lam- 
inar profile  could  be  enough  larger  than  those  for  a turbulent  profile  to  over- 
come the  larger  of  the  entering  air  with  a laminar  profile.  Experimental 

results  for  slots  are  inadequate  to  provide  the  answer  to  this  question.  It 
has,  however,  been  found  experimentally  that  by  proper  shaping  of  the  slot  and 
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diffuser behind  the  slot,  a considerable  portion  of  the  total  be?'!  of  the 
botjtedarr*»layer  air  could  be  preserved , In  any  case,  from  the  above  analysis, 
it. jiay  fyi  seen  that  the  range  of  §§*  for  laminar  and  turbulent  boundary 

i^yerSj’jis  of  the  sane  order  of  magnitude  over  the  middle  range  of  boundary** 
laymr  ^emovfl  ( **  /£)  and  at  least  on  the  basis  of  , no  essential 

, difference  is.  the  mechanism  of  instability  should  be  expected,. 

Static  Stability  of  the  Duct  System 

static  stability  of  a flow  system  is  meant  the  tendency  of  a system 
' when  in  a possible  steady  flow  state  to  return  to  that  state  after  being  slightly 
disturbed  from  it  by  a small  (actually  infinitesimal)  perturbation . The  static 
stability  of  a system  cannot  be  determined  by  a consideration  of  for  the 

slot  alone,  even  when  losses  are  included  in  H.  $ the  characteristics  of  the 
entire  system  must  be  considered.  As  a first  example,  a very  simple  system  con- 
sisting ef  an  entrance  slot  at  a point  of  pressure  R , a diffuser,  a loss-free 
duct  and  an  exit  orifice  also  at  pressure  R . 


The  exit  orifice  can  be  assumed  to  obey  the  simple  pressure  drop  law  (for 
negligible  velocity  in  the  duct) 


u - ^ ys(pi-R)  - - - (9) 

where  G*  is  an  orifice  coefficient.  The  flow  out  the  orifice  qe  is  then 
given  by 

9e-  V2(R  -R)  - - - (10) 


where  K is  a constant  including  both  the  orifice  velocity  coefficient  and 
the  orifice  contraction  coefficient.  The  Intake  characteristics  can  be  given 
by  H,  “ H.-'AH.  , where  H.  is  a function  of  the  boundary-layer  profile  as 
computed  above  and  AH.  is  the  loss  in  total  head  in  the  slot  and  diffuser. 
Roughly,  AH,  may  be  expected  to  be  a function  of  flow  quantity  squared. 

The  flow  in  the  slot  may  be  represented  by  q*  . For  steady  flow,  the  flow 
into  the  system  must  equal  the  flow  out,  and  the  pressures  across  the  system 
must  balance.  Thus 


and 


<le  * % 

H0  - A H,  - R 


1 


— - (11) 
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where  P,  is  total  head  at  the  exit.  This  may  be  written  in  terms  of  qLe  from 
*1.  VO)  as 

p.  - p.  * (l^)*  (18) 

A steady  flow  is  possible  for  some  value  of  <3 » - q*  - <3  when 

H0(<3)  -AH.(q)  = R + <3±  (t&J*  - - - (13) 

Plotting  the  two  sides  of  this  equation  versus  q.  leads  to  a curve  like  that 
shown  below 


It  may  be  seen  that  the  curve  of  R intersects  the  curve  of  H„  — AH.  at  two 
points,  A and  C . Point  A corresponds  to  aero  flow  through  the  system  which 
is  reasonable,  since  there  is  no  static  pressure  difference  across  it  while 
point  C corresponds  to  flow  through  the  system.  Point  B shows  the  even 
larger  flow  through  the  system  which  would  be  possible  if  there  were  no  losses 
at  the  entry,  AH.. 

The  question  now  arises  of  whether  the  flow  which  is  statically  possible 
at  points  A and  G is  also  stable  at  one  or  both  of  these  points.  To  answer 
this  question,  it  is  necessary  to  consider  a small  increase  in  flow  rate,  rSq, 
to  occur.  This  will  cause  an  increase  in  pressure  at  the  entrance  to  the  duct 

of  amount  ^ £ SL  ^ad  811  increase  in  back-pressure  R at  the  exit  of 

amount  . Then  if  H?)  grater  than  , the  mass  of  air 

in  the  duct  will  be  accelerated  toward  the  exit,  thus  further  increasing  the 
flow  through  the  system.  In  this  case  the  steady  flow  is  unstable.  If  on  the 

other  hand,  ■ • I less  than  , the  flow  will  be  decelerated, 

since  the  rise  in  back-pressure  will  be  greater  than  the  rise  in  inlet  pressure. 
In  this  case,  the  flow  is  stable,  since  it  tends  to  return  to  the  original  flow 
rate.  From  the  curves,  it  may  be  seen  that  the  slope,  sKbjarAH.l  is  greater 
than  the  slope,  , for  the  condition  of  no  flow  at  point  A , while  it  is 

smaller  than  at  point  C . Thus,  the  stable  condition  is  that  at  G , with 
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jEIow  through  <,na  system  rathA*  than  that  at  A with  no  flow  though  the  system, 
in  spite  of  th#  fact  that  there  is  no  diffarence  of  static  pressure  across  the 
^Btflmo  It  Is” iapcr'ant  to  note  that  these  conclusions  as  to  stability  depend 
upon  both**the  inlet  characteristics  and  the  system  characteristics,  not  on  the 

inlet  alone?  i.e,*  for  stability,  # It  is  true  that  for 


d(H»  — AH.)  itself  less  than  sero,  the  system  would  necessarily  be  stable  at 
the  sero  flow  condition.  On  the  other  hand,  even  if  positive, 

the  system  can  be  stable  in  the  sero  flow  condition  if  is  sufficiently 

large  there,  which  would  be  the  case  if  the  system  had  large  enough  losses,, 


It  requires  no  particular  change  in  the  above  analysis  to  imagine  that  the 
exit  orifice  is  actually  a second  slot  inlet,  since  for  flow  out  of  the  inlet, 
it  would  have  essentially  the  riov  characteristics  of  an  orifice.  Thus,  the 
interconnection  of  two  identical  slots  by  a duct  might  be  expected  to  lead  to  a 
stable  condition  in  which  flow  enters  one  slot  and  leaves  through  the  other,  the 
condition  of  zero  flow  through  the  system  being  unstable.  The  particular  slot 
rh  acts  in  thb  inlet  would  be  e matter  of  random  choice,  depend?  ug  on  the 
nature  of  initial  disturbances  leading  to  any  given  situation. 


Actually,  the  problem  of  two  identical  slots  discussed  above  is  not  of  as 
much  practical  importance  as  the  one  in  which  there  are  both  two  Inlet  slots 
and  an  exit  orifice  at  reduced  pressure,  which  corresponds  to  the  case  of 
boundary-layer  suction  with  two  slots,  as  shown  belows 


It  should  be  noted  here  that  Norman  and  Holshauser  (ref,  7)  have  analysed 
a similar  system  in  connection  with  twin-duct  intakes  for  engines. 

For  simplicity,  H.(q)  and  AHJq)  are  assumed  to  be  identical  for  the  two 
slots  A and  B , The  flow  rate  through  A ig  <jA  , the  flow  rate  through  B 
is  , and  the  total  flow  is  the  sum  of  and  qB  , The  two  diffusers 
empty  into  a plenum  chamber  in  which  the  flow  velocity  is  assumed  to  be  neg- 
ligibly email.  For  a steady  (but  not  necessarily  stable)  flow  to  be  possible. 
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q - q*  * qB 

H.(q„)  - AHjq*)  — R 
H.(qe)  -AH„(q8)  = R 
R-  Pa  *■  % (*&;) 

These  equations  may  be  represented  graphically  below,  In  the  curves,  due 
account  mist  be  taken  of  the  fact  that  one  or  the  other  of  the  slots  may  act 
as  an  exit  as  well  as  an  intake , lb  may  be  seen  that  the  situation  has  become 


considerably  more  complicated  than  the  one  previously  studied,,  Nevertheless, 
the  seme  basic  principles  apply.  As  in  the  simpler  case  previously  treated,  the 
problem  can  best  be  discussed  in  terms  of  curves  of  H.-AH,.  versus  q and  R 
versus  q , The  exit  orifice  curves  for  three  different  values  of  exit  pressure, 
Pz  , Curves  I,  II  and  III  are  simple  parabolas  corresponding  to  Bq,  (9),  Curve  IV 
for  outflow  through  one  slot  as  an  exit  orifice  is  similarly  plotted  as  a para- 
bola but  for  negative  q , denoting  that  it  subtracts  from  the  total  inlet  fXov0 
Curve  V is  the  basic  inlet  flow  curve  for  a single  slot,  and  has  the  shape  prev- 
iously noted.  Curve  VI  is  the  curve  for  equal  flew  through  the  two  inlet  slots | 
it  is  obtained  by  multiplying  the  flow  quantities  for  Curve  V by  2,  Curves  Vll(a) 
and  VII (b)  bounded  by  points  A,  B ,C  and  Q are  the  curves  for  unequal  and  in- 
flow, and  for  inflow  in  one  slot  and  outflow  in  the  other,  respectively.  Since 
the  plenum  chamber  pressure  is  constant,  the  values  of  H,— AW, , or  P,  (if  the 
flow  ie  out)  must  be  the  same  for  the  two  slots.  Thus,  Curve  711(b)  is  con- 
structed by  subtracting  from  the  inlet  flow  curve  for  one  slot,  V,  the  flow 
quantities  for  the  single  slot  outflow  curve,  IV,  at  the  same  pressure. 

Curve  VH( a)  which  represents  unequal  but  not  reversed  flow  in  the  two  inlets 
is  constructed  by  adding  the  flow  rates  on  the  left  side  of  the  peak  pressure 
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point  5 to  those  on  the  right  side,  giving  segment  A B , Thus,  it  now  app**~* 
that  for  each  of  the  exit  orifice  curves,  there  are  two  solutions  for  plenum 
chamber  pares sure  and  flow  rate,  one  fro*  the  intersection  of  the  exit  orifice 
curves  with  Curve  "VI  (equal  flow)  and  the  other  fro*  an  intersection  with 
Curve  VXI  (unequal  or  reversed  flow).  Of  course,  if  the  exit  pressure  is  low 
enough  to  produce  a total  flow  rate  to  the  right  of  point  A,  there  is  only  one 
possible  flow,  not  two,  corresponding  to  each  inlet  operating  os  Curve  7 to  the 
right  of  point  B, 

It  is  now  necessary  to  inquire  into  the  question  of  which  solution  is  stable 
when  two  exist.  The  criterion  previously  evolved  for  a single  inlet  sysla*  can- 
not be  directly  applied  to  the  curves  of  total  flow  versus  plenum  chamber  pressure 
for  the  entire  system,  and  such  curves,  which,  incidentally,  are  those  directly 
obtained  fro*  experiment,  can  give  no  direct  indication  of  the  stability  of  the 
systs*.  Instead,  it  is  necessary  to  consider  whether  small  disturbances  in  the 
established  flow  rates  in  any  branch  of  the  system  can  upset  the  established 
equilibrium.  If  this  be  done,  it  is  found  analogously  to  the  simple  case  prev- 
iously considered  that  the  conditions  for  stability  are 

3(H,  ap 

3'3a 

. -aa.  <o 

dq_ 

a(H,-AHok  + afHo-AHcL  <Q 

If  the  flow  Is  reversed  in  one  of  the  ducts,  then  the  sign  of  the  derivative 
relating  to  it  in  the  last  equation  is  reversed.  All  derivatives  are  evaluated 
for  the  values  of  the  variables  corresponding  to  the  equilibrium  flow  which  Is 
being  investigated  for  stability  or,  in  other  words,  at  the  values  of  pressure 
and  flow  rate  in  each  branch  of  the  system  for  the  given  case. 

Consider  now  the  case  of  exit  orifice  Curve  II  which  intersects  the  equal 
slot  flow  Curve  VI  at  F and  the  curve  for  reversed  flow  in  one  slot  VXI(b)  at  S, 

In  this  case,  the  point  at  F which  corresponds  to  both  slots  operating  on  the 
single  slot  Inlet  curve  at  point  I is  an  unstable  flow  po?Jit,  since  the  condi- 
tion given  by  2q,  (lit)  is  violated  for  both  slots  having  positive  slopes  of  the 
total  head  versus  flow  rate  curve.  On  the  other  hand,  point  E,  which  corresponds 
to  one  slot  operating  at  point  H on  the  stable  part  of  the  single  slot  curve  and 
the  other  slot  operating  on  point  G of  the  reverssd-ficw  curve,  is  Btable, 

Another  interesting  case  is  the  one  exhibited  by  exit  orifice  Curve  I which 
shows  a stable  point  C at  the  peak  pressure  recovery,  which,  however,  is  obtained 
at  the  expense  of  a siaeable  outflow  fro*  one  slot  operating  at  point  «J  on  the 
outflow  curve.  Thus,  in  general,  the  stable  operating  curve,  if  one  exists,  is 
likely  to  be  of  the  for*  of  Curve  Vli  (points  Q C B A 5)  and  to  involve  .unsym- 
■etrical  flow  and  flow  reversal  in  the  ducts  up  to  point  A-  at  which  point  equal 
flow  in  the  ducts  becomes  the  only  possible  flov0 

Dynamics  of  the  Duct  System 

In  order  to  simplify  the  dynamic  analysis  of  the  air-induction  system, 
consideration  will  be  given  to  motions  which  are  small  perturbations  fro*  a 
steady  flow0  The  system  to  be  considered  consists  of  an  inlet  duct,  a plenum 
chamber  of  volume,  V,  and  an  exit  duct,  as  shown  in  the  figure  below.  Such  a 
system  is  closely  analogous  to  the  Helmhoita  resonator  in  acoustics.  The 
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assumptions  made  for  the  analysis  of  the  Helmhol*  \<k 

resonator  (ref.  8)  nay  be  applied  here.  First,  it  V\\  

is  assumed  that  the  wave  lengths  of  any  oacilla-  \ \ 

tions  which  nay  occur  are  large  compared  to  the  \ \ PLENUM 

dimensions  of  the  system.  Thus,  the  assumption  \ N 

of  unsteady  incompressible  flow  any  be  adopted.  \ o 

Second,  it  is  assumed  that  the  velocities  in  the  \ R — 

plenum  chamber  are  negligibly  snail  and  that  this  \ 

chamber  nay  be  considered  to  be  an  air  spring.  The  j 

inertia  effects  in  the  system  are  assumed  to  be  1 

derived  from  the  masses  of  air  in  the  inlet  and  exit 

ducts  and  from  any  apparent  mass  contributions  due  to  the  air  at  the  esds  of  these 
ducts.  The  volume  flow  rates  are  q,  and  qe  for  the  inlet  and  outlet  ducts  respec- 
tively. The  total  pressure  at  the  entrance  to  the  inlet  duct  is  denoted  by  Hs  j 
thus,  total  pressure  is  assumed  to  include  the  effect  of  all  losses  in  the  duct 
itself  as  well  as  the  boundary-layer  and  entrance  losses.  The  pressure  immediately 
ahead  of  the  exit  duct  is  equal  to  the  plenum  chamber  pressure  Pc . 

The  spring  constant  of  the  plenum  chamber  can  now  be  computed.  The  mass  incre- 
JM*it  accumulated  in  the  plenum  chamber  due  to  incremental  flow  rate  changes  q*  and 
q,  is 


SMc  - 


(15) 


The  corresponding  change  in  density  is 


6^  - = % \ 


(16) 


where  Y is  the  volume  of  the  plenum  chamber.  The  pressure  is  related  to  the 
density  by 


n'a  s 4B 

a Jp 


where  a is  the  * velocity  of  sound.  Then  rff?  , or 

fiR  - a1  % J (q,  - qe ) d* 
The  tern  a*(J^  will  be  called  It  , so  that 


6R  - 


(17) 


- - - (18) 


(19) 


The  flow  in  the  inlet  duct  will  next  be  considered  as  a one-dimensional,  un- 
steady, incompressible  flow.  The  Euler  equation  is 


— • * (flj)' 


Integrating  this  over  the  length  of  the  duet,  .£  , gives 


j P §£  ^ -o  (21) 

*'•  •/»  Ja 
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The  second  integrand  is  an  exact  differential,  3%  that  the  last  equation  may 
be  written 


O 


or 


dx 


+ R - Hs  - O 


- “ (22) 


since  ( (*£■*+  p)  is  Hs  at  tho  duct  entrance  and  R.  at  the  duct  exito  Since 
the  flow  is  incompressible,  the  velocity  along  the  inlet  duct  is  given  by 


u = qyA 

where  A is  the  cross-sectional  area  of  the  duct  at  any  point „ Thus,  the 
Eulerian  equation  may  be  written  as 


? If*  + R'Ha=0  - - - (23) 

Using  As  as  the  area  of  the  slot  entrance,  this  may  be  nondimensionalized 
to  give 


§?*  P^* t + PJ-K-O 

The  remaining  integral  is  a dimensionless  function  of  duct  geometry  only# 
Using 

■%  = f .d..feil 

4,  (%s  ) 

the  equation  of  motion  is 


<&) 

(25) 


H,  * O - » - (26) 

For  the  unperturbed  motion  Ro  — ™ O , It  iB  therefore  only  necessary 

to  consider  the  changes  6 R and  <5  Ha  arising  from  the  perturbation  flov0 
<5R  is  given  by  Eq„  (19),  while  for  small  perturbations,  6 Hs  is  given  by 


fN»-  %% 


since,  as  before 


■>,  5 _ -presents  only  the  perturbation  flown 

slope  of  the  characteristic  curve  for  the  inlet  slot,  which  has  previously 


(27) 

is  ths 


d C - A Ho)  ] 


been  discussed  in  connection  with  static  instability,  [-^q4 

Combining  Eqs.  (19),  (26)  and  (27)  leads  to  the  following  equation  for 
the  inlet  duct  flow* 


- (28) 
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For  the  outlet  duct*  ft  very  similar  aquation  can  be  obtained,  sanely. 


where  JU  3 % , and  ,Aft  are  now  to  be  computed  for  the  outlet  duct  and 

is  the  slope  *of  the  characteristic  curve  of  the  outlet  duct  which  has  been 
also  discussed  previously  in  the  section  on  static  instability*  Now  let 


(29) 


and 


?•-  1*1* dt 

?.  - 


Eqs.  (28)  and  (29)  then  become 

-o 


(30) 


(31) 


These  may  be  recognized  as  analogous  to  the  equations  of  motion  for  two  masses 
covqpled  by  a spring  and  having  dampers  attached*  A positive  value  of  ^4* 
may  be  seen  to  correspond  to  negative  damping  of  one  of  the  masses  while  ^ 

a positive  value  of  corresponds  to  positive  damping  of  the  other  mass* 

The  system  is  illustrated  below. 


MASS 

MASS 

II 

SPRING 

n 

h—j  l-q 


NEGATIVE  DAMPER 


POb/TIVE  DAMPER 


Helmholtz  Resonator  Considerations 


If  there  is  no  initial  flow  through  the  system  so  that  and 

are  effectively  zero,  Hqe.  (31)  reduce  to  the  equations  for  a Helmholtz  reson- 
ator with  two  nadks  (ref.  8 )•  The  equations  may  then  be  written  as 


-r  uV^9  -*),*§*  *°  1 


(32) 
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where 


and 


“•  - ^fr. 


- ai 


(33) 


The  formulas  for  to*  and  ^ are  the  usual  formulas  for  simple  Helmholtz 
resonators  having  one  neck  (ref » 8),  except  for  the  factors  and  -^e 
which  are  corrections  for  the  shape  of  the  ducts  or  necks  in  the  resonator 
analogy.  Eqs.  (32)  may  be  solved  to  give  the  frequency  of  the  resonator 
with  two  necks  in  terms  of  the  frequencies  it  would  have  with  each  neck  alone. 
This  frequency  (*>  is  given  by  the  solution  of 


which  leads  to 


-60$  (-U>z*ag) 


= O 


(3U) 


tO4  -u)*(o)£  + 6t)£)  “O 


(35) 


This  equation  has  the  roots 


to1  «*•  o 


and 


03"*=  tA?  oO£ 


- - - (36) 


The  root  a)  - o corresponds  to  the  fact  that  air  can  leak  slowly  into  one 
neck  and  out  the  other,  and  is  of  no  interest  in  the  present  problem.  The 
frequency  of  interest  is  the  one  corresponding  to  the  second  root 

■u)  - + cot  (37/ 


or 

*’  + A.)  ---(38) 

\ 

If  the  exit  duct  length,  ,£e  , is  large  compared  to  the  inlet  duct  length,  £a  , 
the  term  is  small  compared  to  > and  frequency  is 

approximately  equal  to  the  frequency  of  a resonator  having  only  the  inlet  duct 
as  a neck.  This  was  the  case  in  the  experiments,  which  will  be  discussed  later 
in  this  report. 

For  calculations  which  are  made  to  afford  a comparison  of  theoretical  and 
experimental  results,  it  is  necessary  to  compute  the  value  of  -€,s  for  the  duct 
shape  of  interest.  In  the  present  experiments,  the  duct  was  a two-dimensional 
diffuser  with  a linear  expansion  of  area  from  the  slot  to  the  plenum  chamber. 
Taking  the  ratio  of  areas  at  the  two  ends  to  be  n , the  area  at  any  section 
a distance  x from  the  slot  is  given  by 

A-As[lt(n-,)J]  (39) 


WADG  TR-53-189 


- 13  - 


The  integral  for  §r  ia  then  given  by 

*•“  [ life 

Substituting  this  into  Eq.  (33)  gives  for  the  frequency 


I n (r. ) 
n - 1 


- - - (Uo) 


to 


M A,  (n-l) 

s - a 


— (Ill) 


Dynamic  Stability  Conditions 


The  linearized  differential  equations  of  motion  for  small  perturbations 
of  the  duct  system,  Eqs.  (31)  way  be  analyzed  to  determine  the  criteria  for 
dynamic  stability  of  the  system.  Since  the  stability  is  defined  in  terms  of 
response  to  infinitesimal  disturbances,  these  stability  criteria  based  on  the 
small-perturbation  equations  are  exact  to  the  extent  that  the  physical  para- 
meters involved  in  the  coefficients  of  the  differential  equation  can  be  deter- 
mined exactly.  The  solutions  to  the  differential  equations  which  describe  the 
motion  of  the  system  are,  on  the  other  hand,  not  exact  for  non-linear  systems, 
such  as  those  being  considered  here.  The  accuracy  of  the  calculated  frequency, 
for  instance,  depends  on  the  amplitude  of  any  motions  described  being  small 
enough. 

The  equations  of  motion,  Eqs.  (31)*  may  be  written  as 


4 bgd~>e  -t-  (Of  * O 


(hs) 


where 


bK 


b« 


- - 1£  (pfa) 


? 


m. 

irw, 


(i>3) 


The  aolutions  are  of  the  form 


F-t 


“ c.e 

?«  - c,e-  j 

and  the  characteristic  equation  is  then 

(p‘^P*4P!*b«P^i  - -O 

or 


— m 


— - (h$) 


P4  4 ( h*  * be  )pa  + (io‘  u>*  te)p0  + (coi  be  +•  av2  bjJ  P ~ O 
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Factoring  out  the  root  P»o  leads  to  the  cubic 

/ . 

P * * (b^fajp*  + (cO*  +'J>2  + b*  ba)P  + (o)?b,  * cq?b»)  -O  (1*6) 


According  to  Routh's  criteria  (ref.  9)*  P will  not  hare  a positive  real  part 
(i.e.,  the  system  will  be  stable)  if  all  the  coefficients  of  the  equation  are 
positive  and  if,  for  the  equation 


P*  + BP*  ♦ Cp  D — O t BC  - D > O (Routh‘8  Discriminant) (1*7) 

or  in  this  case 

(bs  * tfe  + b*)  - ( oDa  bs ) > O 

This  can  be  simplified  to 

CO*  b,  + tOe  + - - - (3*8) 

The  requirement  that  the  coefficients  of  the  characteristic  equation  all 
be  positive  is  sufficient  to  insure  static  stability  of  the  type  discussed 
previously.  In  particular,  the  requirement  that  the  constant  coefficient  of 
the  equation  must  be  positive  leads  to 

cO*  * cog  bs  ~^0 

or 

[~fc  If.  “■  >0  — Cte) 

This  finally  is  reduced  to 

~ + f|e  ^ ° - - - (50) 

which,  with  due  account  of  changes  in  notation,  is  equivalent  to  Eqs.  (ll*)  for 
static  stability  in  the  same  case.  Since  a real  cubic  has  either  three  real 
roots  or  a real  root  and  two  complex  conjugate  roots,  and  since  the  constant 
coefficient  is  equal  to  minus  the  product  of  the  roots,  a negative  value  of 
this  coefficient  indicates  the  existence  of  at  least  one  real  positive  root. 

This  represents  a non-oscillating  divergent  type  of  instability  which  is  what 
would  be  expected  in  the  case  of  static  instability. 

If  the  coefficient  of  P*  is  negative  but  the  constant  coefficient  is  not, 
(i.e.g  the  system  is  statically  stable)  the  instability  is  oscillatory  diver- 
gent. This  is  clear  from  the  fact  that  the  coefficient  of  P*  is  equal  to 
the  sum  of  the  roots  of  the  cubic,  so  that  if  there  is  no  positive  real  root, 
the  coefficient  can  only  be  negative  if  the  real  part  of  the  pair  of  eoaplax 
conjugate  roots  is  positive.  The  requirement  is 

b4  * be  >-0 


or 
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> o 


(51) 


i 


i 


This  is  similar  to  the  condition  for  static  stability  except  that  the  damping 
term  for  each  duct  is  divided  by  the  mass  term  for  that  duct.  This  is  connected 
with  the  fact  that  in  the  analogous  oscillatory  motion  of  two  masses  coupled  by 
a spring,  the  amplitudes  of  the  two  elements  are  inversely  proportional  to  their 
masses.  Thus,  a damper  on  the  larger  mass  is  less  effective  than  one  on  the 
smaller  mass.  Similarly,  in  the  present  case,  high  exit  losses  ( §£  large), 
which  may  be  sufficient  to  provide  static  stability  according  to  Eq.  (00)  may  be 
ineffective  in  providing  dynamic  stability  according  to  Eq.  (01),  if  the  effec- 
tive exit  mass,  , is  large  compared  to  the  effective  inlet  mass. 

Finally,  if  all  the  coefficients  of  the  characteristic  equation  are  made 
positive,  the  sole  criterion  for  stability  is  Routh’s  discriminant,  Eq.  (U8). 
Substituting  the  values  of  ^ , <4.  , bs  and  be  into  the  first  two  terms  of 
this  equation,  there  results 


3P, 


6u) 


+ b5  b<-  ( b6  + be)  >•  o 


(52) 


Thus,  if  is  positive  ( bs  negative)  while  — 

(static  stability  condition)  and  (bs  + bt ) is  also  positive,  this  condition 
can  only  be  satisfied  if  the  first  term  is  sufficiently  positive.  In  order 
for  this  term  to  be  positive  at  all,  the  sum  of  the  damping  coefficients  divided 
by  the  squares  of  the  mass  coefficients  must  be  positive.  This  is  an  even  more 

severe  requirement  for  large  (large  exit  loss)  than  Eq.  (01),  if  the 

effective  mass  of  the  exit  duct  is  relatively  large.  If  the  term  in  brackets 
is  positive,  then  increasing  k (which  amounts  to  decreasing  the  plenum  volume) 
is  favorable  for  dynamic  stability. 

Non-Linear  Effects 


From  the  stability  criteria  given  above,  it  may  be  seen  that  if  the  effec- 
tive mass  of  the  exit  duct  is  large,  a statically  stable  system  can  easily  be 
dynamically  unstable.  According  to  the  linearized  equations  given  above,  the 
amplitude  of  oscillation  would  increase  in  time  without  limit.  Actually,  due 
to  the  nonlinearity  of  the  system,  the  oscillations  will  settle  down  to  a 
periodic  motion  at  some  finite  amplitude.  This  motion  is  called  a "limit  cycle" 
in  nonlinear  mechanics.  This  is,  of  course,  what  is  actually  observed  in  prac- 
tical cases  of  dynamic  instability.  To  investigate  qualitatively  the  nature  of 
the  limit  cycle,  it  is  convenient  to  consider  a case  in  which  the  effective  mass 
of  the  exit  duct  is  so  high  that  it  experiences  relatively  little  amplitude  of 
oscillation  compared  to  the  inlet  duct.  In  this  case,  the  motion  may  be  con- 
sidered as  that  of  the  inlet  duct  only.  The  equation  of  motion  is  then 


where  bs(5«,  <2.) 
taneous  value  of 
be  multiplied  by 
results 


, q0)%„  -hcO*€,5=0  - - - (. 

is  a function  of  the  mean  flow  rate  and  the  Instan- 
ts, and  is  considered  to  be  constant.  If  this  equation 
and  integrated  over  a period  of  the  motion,  T,  there 


+ 


♦ cDa 


(0U) 
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This  is  equivalent  to 


or 


js.4.  - {%*)*& 

f Hif + <**  ]o  ♦ ^ [it , «j.)  (f.fdt ' - o 


The  first  tern  is  zero  if  T is  the  period,  so  the  condition  for  periodic  notion 
is 

j^k>s(5*>  3®)(§s)  d+  =o  - - - (55) 

I 

This  is  what  night  have  been  expected;  it  states  that  the  average  energy  input 
to  the  systen  over  a cycle  must  be  zero.  The  condition  nay  be  written  as 

J -o  ---(56) 

If  bs  is  small  compared  to  to5  , the  solution  of  Bq.  (53)  will  be  very 
nearly  sinusoidal,  namely. 


“ A S;«  cD3t 


and 


Acdf  Cos 


The  curve  of  ^5^5,  versus 
the  figure  below. 


The  integral  will  be  equal  to 
-jjf*  multiplied  by  the  value 


will  then  be  a circle  of  radius  A as  shown  in 


the  limit  of  the  sum  of  each  elementary  strip 
of  ii)j  bj  appropriate  to  that  value  of  . 


As  a very  simple  case,  consider  a curve  of  of  the  type  shown  below. 

The  slope  is  a positive  constant  to  the  left  of  Ck  and  a negative  constant 
to  the  right. 
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Here  q6  is  the  wear,  flow  rate  and  qc  is  the  point  on  the  curve  where 
changes  from  positive  to  negative*  Then  at  s ^changes  from 

positive  to  negative o Thus,  on  a diagram  of  %^,/nX,  versus  %%  s the  limit 
cycle  corresponds  to  the  case  when 

bSiA,  * bSfiAa  *0  - - - (57) 

0 s • 4 

where  b&)  is  the  value  of  b,  for  %,<  %ft<.and  is  the  value  for  > f>Sc 
while  A,  and  Aa  are  the  areas  shown  on  the  diagram  below0 


It  is  obvious  from  these  considerations  that  in  some  cases  oscillation  is 
possible  even  when  qQ  is  large  enough  so  that  is  negative.  In  this 

case,  the  flow  would  be  dynamically  stable  for  small  disturbances,  but  a large 
disturbance  of  the  right  kind  could  start  an  oscillation  of  finite  amplitude. 
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DESCRIPTION  OF  MODELS  AND  EXPERIMENTAL  APPARATUS 


The  experimental  work  was  carried  out  in  the  Cornell  Aeronautical  Laboratory 
3*  by  16*  subsonic  induction  tunnel,  which  is  of  the  non-return  type.  This  tunnel 
is  powered  by  the  600  horsepower  Nash  poops  used  for  pressurlsation  of  the  Cd£eL. 
££'  by  12*  Variable  Density  Tunnel  which  are  used  to  supply  the  elector  jets  for 
thic  tunnel.  Most  of  the  tests  were  mads  using  the  turbulent  boundary  layer  on 
the  16*  wall  of  this  tunnel  as  the  induction  air0  A slot  13o?5*  wide  was  wade 
in  the  wall  of  the  tunnel  and  opened  into  a diffuser  section  leading  to  the 
plenum  chamber  . The  boundary-layer  air  was  evacuated  from  the  plenum  chamber 
through  a line  connected  to  a 1*0  horsepower  Kinney  positive  displacement  vacuum 
pusqpo  A bypass  valve  was  also  provided  in  this  line  to  allow  auction  of  air  from 
the  room  along  with  that  fro*  the  plenum  chamber 0 This  bypass  permitted  a wider 
range  of  suction  pressures  and  volume  flows  than  was  possible  with  a direct  con- 
nection to  the  pump0  It  also  permitted  dynamic  conditions  at  the  exit  orifice 
which  were  not  as  strongly  dependent  on  the  vacuum  pump  characteristics  as  with 
a closed  line.  The  general  arrangement  of  the  wind  tunnel,  plenum  chamber,  and 
vacuum  line  is  shown  in  Figs.  1 and  2. 

The  slot  entry  and  diffuser  were  designed  to  provide  good  pressure  recovery 
of  the  entering  boundary-layer  air0  Based  on  tests  by  PLerpont  (ref.  10),  the 
slot  shape  and  diffuser  design  shown  in  Fig.  3 were  selected.  The  area  ratio  of 
the  diffuser  between  plenum  chamber  entrance  and  throat  was,  however,  much  larger 
than  that  used  in  PLerpont f a tests  in  order  to  avoid  the  loss  of  dynamic  head 
which  would  otherwise  occur  at  the  plenum  entrance.  The  diffuser  was  two- 
dimensional  and  was  formed  by  two  blocks.  The  forward  block  was  fixed,  while 
the  aft  block  was  moveable  to  provide  for  variable  slot  width.  The  diffuser 
area  ratio  varied  fro*  16  at  a elot  width  of  3/16"  to  7.5  at  a slot  width  of 
7/16".  By  contrast,  PLerpont 's  area  ratio  was  2.0.  The  diffuser  half-angle 
remained  fixed  at  6 , since  the  displacement  of  the  aft  block  involved  no  rota- 
tion ~f  the  wall. 

Total  head  surveys  were  made  to  obtain  the  boundary-layer  velocity 
profiles  2.75"  ahead  of  the  slot  along  the  wall  centerline.  Also,  total-head 
tubes  and  re  verse- total-head  tubes  were  located  in  the  diffuser  at  four  stations 
along  the  slot  span  to  check  flow  uniformity  in  the  slot.  This  is  shown  in 
Fig.  lu  Tunnel  static  pressures  were  measured  by  means  of  wall  static  taps, 
and  plenum  chamber  pressure  was  measured  at  several  points  in  the  plenum 
chamber.  These  pressure  readings  were  indicated  on  a manometer  board  and 
recorded  photographically.  The  indications  were  relatively  steady  in  the 
presence  of  high-frequency  oscillations  because  of  the  damping  effect  of  the 
long  pressure  lines  and  the  manometer  itself  <,  Mass  flow  from  the  plenum 
chamber  was  measured  by  means  of  a calibrated  Venturi  tube. 

Oscillating  pressures  ware  measured  by  means  of  a Stathas  Type  P-6  dynamic 
pressure  pickup  of  the  strain  gage  type  and  recorded  on  a Brush  oscillograph. 

•”  various  times,  the  dynamic  pickup  v^e  vented  to "various  pressure  taps  in  the 
system  in  order  obtain  a sampling  of  dynamic  data  at  all  points.  Near  the 
end  of  the  program  the  Brush  oscillograph  was  unavailable  and  some  data  were 
taken  with  a carbon  microphone,  anplifier  and  cathode-ray  oscilloscope.  This 
arrangement  turned  out  to  be  less  satisfactory  because  of  its  sensitivity  to 
noise  in  the  tunnel  and  compressor  system.  All  dynamic  data  presented  in  this 
report  is  based  on  measurements  with  the  pressure  pickup. 
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la  order  to  experiment  with  laminar  boundary- layers  * a second  model  was 
constructed.  This  was  of  airfoil  shape  with  a flat  lover  surface.  The  suction 
slot  and  diffuser  were  scaled-down  versions  of  those  used  in  the  wall  model, 
and  were  located  12  inches  from  the  airfoil  leading  edge.  The  toads 1 spanned 
the  tunnel  along  the  16-inch  dimension  and  had  its  lowex-  surface  3/U  inch 
from  one  -of  the  walls.  The  aft  portion  of  the  model  was  used  as  the  plenum 
chamber.  The  plenum  chamber  was  connected  to  the  gate  valve,  Venturi  tube  and 
vacuum  pump  through  a pair  of  3/U  inch  I.D.  rubber  tubes  connected  at  each  side 
of  the  airfoil.  A schematic  diagram  of  the  test  set-up  is  shown  in  Fig.  £«, 
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EXPERIMENTAL  RESULTS  AND  DISCUSSION 


Tunnel  Wall  Boundary  Layer-Wide  Span  Slots 

The  first  tests  conducted  in  this  program  were  made  with  13.75  inch  span 
slots  on  the  tunnel  wall.  Typical  houndary-layer  velocity  profiles  ahead  of 
the  slot  are  shown  in  Fig.  7.  Both  static  and  dynamic  instability  were  en- 
countered in  these  tests  for  flow  rates  lower  than  that  required  to  suck  the 
entire  boundary  layer  across  the  slot  span.  The  static  instability  took  the 
form  of  a nonuni form  distribution  of  flow  across  the  diffuser.  At  low  flow 
rates,  the  flow  actually  reversed  at  one  end  of  the  slot,  A typical  set  of 
results  is  shown  in  Fig.  8,  In  this  case,  the  boundary  layer  thickness  was 
about  .68  inches.  The  results  shown  are  for  suction  flow  rates  corresponding 
to  sucking  .Oh  inches  and  .2  inches  of  the  boundary  layer.  In  a crude  way,  the 
static  stability  analysis  given  earlier  for  a two-slot  system  may  be  applied  to 
a system  with  a single  wide  slot  by  assuming  it  to  be  divided  in  two.  It  is 
then  clear  that  the  more  stable  flow  at  low  flow  rates  would,  in  fact,  be  one 
in  which  the  flow  rates  in  the  two  sides  of  the  system  is  unequal.  The  asymmetry 
of  the  flow  along  the  span  decreased  as  the  suction  flow  quantity  increased. 

This  effect  is  better  illustrated  for  a typical  case  in  Fig.  9 in  terms  of  the 
asymmetry  index  J defined  as 


t — Ln 

J ~ 

where  umw  and  um£„  are  the  maximum  and  minimum  velocities  measured  across 
the  diffuser,  Q is  the  total  volume  flow  rate  through  the  system,  and  A is 
the  area  of  the  diffuser  at  the  section  where  the  measurements  were  made.  It 
may  be  seen  that  for  flow  rates  approaching  half  those  necessary  to  suck  the 
entire  boundary  layer,  the  asymmetry  becomes  small.  According  to  the  theory  for 
two-slot  systems  previously  presented,  this  should  occur  when  the  peak  of  the 
plenum  pressure  recovery  curve  has  been  passed  for  the  flew  in  both  parts  of  the 
slot.  It  proved  to  be  difficult  to  verify  this  theoretical  requirement  because 
the  details  of  the  three-dimensional  flow  In  the  slot  do  not  fit  the  simplifying 
assumptions  of  the  two-slot  theory,  and  because  whenever  this  type  of  static  in- 
stability cc cured  in  the  present  tests,  it  was  also  accompanied  by  dynamic  instab- 
ility with  attendant  oscillations.  It  was,  however,  observed  that  the  diminution 
of  asymmetries  occurred  at  flow  rates  where  the  slope  of  the  curve  of  plenum  chamber 
pressure  versus  flow  rate  was  negative. 

As  mentioned  above,  the  wide-span  slots  all  showed  both  static  instability 
and  dynamic  instability  simultaneously,  A typical  plenum  chamber  pressure  record 
made  with  a dynamic  pressure  pickup  is  shown  in  Fig.  11 j it  may  be  seen  that, 
although  the  record  is  fairly  regular  and  periodic,  it  is  not  sinusoidal.  This 
is  evidence  of  the  nonlinearities  in  the  system  and  was  to  have  been  expected. 

In  spite  of  this  fact,  the  Helmholtz  resonator  frequencies  calculated  on  the 
basis  of  the  linearized  theory  (Eq.  Ul)  agreed  fairly  well  with  measured  fre- 
quencies in  most  cases,  A comparison  of  calculated  and  measured  frequencies  is 
given  in  Table  I, 
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TABLE  X 

COMPARISON  OF  CALCULATED  AND  MEASURED  OSCILLATION  FREQUENCIES 


SLOT  SPAN 

SLOT  WIDTH 

CALCULATED 
FREQUENCY  (cps) 

MEASURED 
FREQUENCY  (cps) 

13.75  ins. 

3/16  in. 

30 

27.0  i.5 

13.75  ins. 

5/16  in. 

33 

28.5  i.5 

13.75  ins. 

7/16  in. 

37 

30.0  -.7 

7.85  ins. 

5/16  in. 

25 

21.5  i.5 

3.85  ins. 

5/16  in. 

18.5 

13  -20 

The  calculated  frequencies  according  to  Eq.  (ip.)  are  based  on  the  assumption 
that  the  effective  mass  of  the  exit  duct  is  very  large  compared  to  the  effective 
mass  of  the  inlet  duct.  If  this  were  true,  the  frequency  of  the  system  with  a 
7*85  in.  span  slot  should  b ( * 1.32  ) times  the  frequency  of  the  same  system 
with  a 13.75  in.  span  slot.  The  ratio  of  the  measured  frequencies  from  Table  I is 
^^§—1 32,  which  provides  a good  check  of  this  assumption.  The  uncertainty  of  the 
measured  frequency  is  associated  with  a generally  small  rise  in  frequency  with  flow 
rate.  This  effect  became  most  pronounced  for  the  3*85  in.  span  slot  listed  in 
Table  I for  which  the  frequency  varied  from  13  to  20  cps.  This  effect  is  another 
indication  of  the  nonlinearity  of  the  system.  In  order  to  check  on  the  possibility 
that  the  flow  in  the  wind  tunnel  might  be  coupled  to  the  plenum  chamber  oscillation, 
teste  were  run  with  varying  lengths  of  tunnel,  but  the  effects  were  generally  small. 

The  ranges  of  dynamic  stability  observed  in  the  tests  of  the  wide-span  slots 
are  shown  in  Fig.  10  on  the  curves  of  plenum  chamber  pressure  versus  flow  rate. 

Test  points  on  these  curves  are  shaded  if  oscillation  occurred.  Also  shown  on 
this  figure  is  a curve  of  the  boundary-layer  profile  and  a curve  of  flow  rate 
through  the  system  versus  amount  of  boundary  layer  which  would  be  removed  if  the 
flow  were  uniform, in  terms  of  distance  from  the  wall  of  the  outermost  streamline 
entering  the  slot.  It  may  be  observed  that  for  all  slot  widths,  the  plenum  pres- 
sure drops  with  flow  rate  for  small  flow  rates  and  the  system  is  at  first  stable. 

The  reasons  for  this  effect  have  not  been  determined.  It  may  be  that  for  small 
flow  rates  with  large  asymmetry  along  the  slot, mixing  losses  in  the  diffuser  and 
separation  of  the  low-energy  air  from  the  diffuser  wall  are  responsible.  On  the 
other  hand,  it  may  be  that  this  effect  iB  simply  a dynamical  one,  due  to  the  fact 
that  with  asymmetrical  flow  such  as  might  be  encountered  in  this  flow  regime,  the 
losses  associated  with  the  reversed  flow  are  sufficient  to  stabilize  the  system. 

The  data  at  hand  are  insufficient  to  settle  this  question.  In  any  case,  for  all 
slot  widths,  the  oscillations  do  not  start  until  the  slope  of  the  curve  of  plenum 
chamber  pressure  versus  mass  flow  has  become  positive,.  The  oscillations  do  not 
necessarily  stop,  however,  when  the  slope  has  become  negative.  This  should  not 
be  surprising,  since  as  on  the  segment  BC  of  the  curve  for  the  two-slot  model. 


W ADC  TR-53-189 


discussed  cn  page  8 , the  plenum  Chamber  total  pressure  may  shew  a falling  pres- 
sure characteristic  while  one  or  the  other  of  the  slots  is  operating  on  a rising 
pressure  characteristic.  If  any  element  of  the  system  has  such  a positive  slope 
of  the  total  head  versus  mass  flow,  then  instability  is  possible.  Since  a wide- 
span  single  slot  is  potentially  a multiple-slot  system  if  s panwise  variations  of 
flow  rate  occur,  the  two-slot  model  may  provide  some  insight  into  its  operation. 
In  fact,  ignoring  the  initial  pressure  drop  in  the  curves,  the  curves  of  Pig.  ID 
show  a certain  vague  similarity  to  the  theoretical  two-slot  svstem  curve  on 
page  8 , in  that  the  plenum  pressures  show  two  maximum  points.  The  differences 
between  the  two-slot  model  and  a wide-span  slot  are  not  negligible,  since  mixing 
between  the  air  entering  and  leaving  the  system  can  take  place  at  high  velocity 
in  the  diffuser  instead  of  at  negligible  velocity  in  the  plenum  chamber.  Thus, 
the  plenum  chamber  pressure  should  be  lower  than  would  be  expected  on  the  basis 
of  the  total  head  of  the  entering  air  and  estimated  slot  losses.  Further,  even 
if  the  slot  flow  were  truly  two-dimensional,  it  would  not  be  known  how  the  losses 
in  unsteady  flow  are  related  to  the  steady-flow  losses. 

In  Fig.  12  are  presented  curves  of  plenum  pressure  versus  flow  rate  for  a 
3/l6  in.  wide  wide-span 1 slot  with  a throttling  device  in  the  diffuser  (as  shown 
in  Fig.  5) ♦ Here  again,  points  at  which  instability  was  encountered  are  shaded. 

It  may  be  seen  that  ranges  of  dynamic  instability  occur  whenever  the  curve  of 
plenum  chamber  pressure  versus  flow  rate  has  a positive  slope.  With  the  throt- 
tling valve  full  closed  (which  still  allowed  passage  of  some  flow),  this  slope 
was  always  negative  and  no  instability  occurred. 

Tunnel  Wall  Boundary  Layer  - Slots  with  Splitters 

In  an  attempt  to  improve  the  flow  uniformity  and  stability  of  the  wide-span 
slots,  the  slot  entry  and  diffuser  were  divided  into  four  equal  sections,  each  of 
3.85  in.  span,  by  means  of  splitter  plates.  The  plenum  was  not  divided,  however. 
Tests  of  this  arrangement  showed  no  dynamic  instability,  but  there  was  static  in- 
stability? evidenced  by  nonuni form  flow  rates  through  the  various  sections  of  the 
duct.  In  fact,  for  very  small  flow  rates,  all  of  the  flow  tended  to  go  through 
one  section  of  the  duct  with  no  flow  or  a small  reverse  flow  through  the  other 
sections.  Various  sections  of  the  duct  would  be  operative  for  one  total  flow 
rate  and  inoperative  for  other  flow  rates  in  a random  fashion.  Once  the  flow  rate 
was  equal  to  the  flow  necessary  to  suck  the  entire  boundary  layer  across  the  spaa 
of  the  duct,  all  sections  became  operative  and  uniform  flow  was  approached.  Typical 
results  for  measured  diffuser  velocities  in  the  four  sections  jf  the  duct  are  riven 
in  Table  .II. 
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TABLE  II 


DIFFUSER  VELOCITIES  FOR  VARIOUS  FLOW  RATES  - FT. /SEC. 
Tunnel  V«  270  ft/sec  Q ~ ft3/sec 


Q 

0.857 

2.00 

3.78 

6.05 

A 

15U  1 hh  163  159 

155  20li  201  206 

180  210  202  161 

185  170  155  126 

n 

0 0 0 0 

0 0 0 0 

112  150  11*6  158 

150  lh9  178  166 

c 

0 0 0 0 

151  162  167  l6l 

0 0 0 0 

165  lii5  180  l6o 

D 

0 0 0 0 

0 0 0 0 

litl  175  189  206 

118  175  0 196 

It  may  be  noted  that  a quantity  of  flow  in  an  active  section  tended  to  be  approxi- 
mately equal  to  the  flow  which  would  just  take  in  all  of  the  boundary-layer  air 
ahead  of  the  section  (one-fourth  of  the  total  boundary  layer  air  ahead  of  the 
entire  slot).  Thus,  at  Q - .857  (ft^/sec),  one  section  was  active;  at  Q = 2.00, 
two  sections  were  active;  at  Q - 3.78,  three  sections  were  active;  and  at  Q » 6.05, 
which  is  just  about  enough  to  remove  the  boundary-layer  air  ahead  of  the  slot,  four 
sections  were  active. 


It  was  desired  to  verify  that  the  elimination  of  the  dynamic  instability  was 
not  simply  due  to  the  increased  losses  arising  from  skin  friction  on  the  splitter 
plates,  although  this  seemed  unlikely.  This  could  easily  be  done  by  sealing  three 
of  the  four  sections  of  the  duct  and  applying  suction  to  the  fourth.  It  was  found 
that,  with  three  sections  sealed,  dynamic  instability  reappeared  with  oscillations 
at  a frequency  appropriate  to  the  reduced  slot  area.  Curves  of  plenum  chamber  pres- 
sure versus  flow  rate  are  presented  in  Fig.  13  for  a 5/16  in.  slot  with  the  splitters, 
without  splitters,  and  with  splitters  and  three  sections  sealed.  The  curves  clearly 
show  that  effect  of  the  splitters  on  dynamic  stability  is  not  an  effect  due  to  in- 
creasing duct  losses,  but  is  a multiple-slot  effect.  Apparently  what  happens  is 
that  the  flow  is  statically  stable  only  under  assymetrical  flow  conditions  which 
call  for  almost  complete  boundary-layer  suction  through  operative  sections  of  the 
duct.  Under  these  conditions,  the  flow  is  dynamically  stable. 

Airfoil  Model  Boundary  Layer 

In  an  effort  to  obtain  experimental  results  on  the  stability  of  a duct  system 
ingesting  a laminar  boundary  layer,  tests  were  run  with  the  airfoil  model  described 
above  (Fig.  6)  at  tunnel  speeds  between  iiO  and  100  ft/sec  both  with  and  without 
stilling  screens  in  the  tunnel.  A slot  width  of  ,05  in.  was  employed  correspond- 
ing to  the  anticipated  boundary-layer  thickness. at  100  ft/sec.  It  was  found  that 
neither  static  nor  dynamic  instability  could  be  observed.  Tests  were  then  also 
made  at  100  ft/aec  with  a transition  wire  on  the  airfoil  to  produce  a turbulent 
boundary  layer  and  again  no  instability  was  observed.  It  appears  that  this  stab- 
ility, as  compared  with  the  larger  model,  may  have  been  due  to  (1)  high  slot  and 
duct  losses  at  the  small  scale  of  these  tests,  (2)  relatively  high  exit  losses  as 
compared  to  the  larger  model  or  (3)  relatively  low  plenum  chamber  volume  compared 
to  the  larger  model.  Within  the  limited  "cope  of  the  present  program,  it  wac 
possible  to  embark  on  the  comprehensive  tests  necessary  to  determine  the  effects  of 
all  these  variables.  Items  (2)  and  (3)  could  be  checked  within  a limited  range  by 
using  a partial-span  slot  to  reduce  the  inflow  to  the  system.  The  slot  span  was 
reduced  to  one-fourth  of  ite  original  li;  in.  value,  but  the  system  remained  stable, 
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Theoretical  and  experimental  studies  of  elementary  air-induction  systems 
for  boundary-layer  auction  lead  to  the  following  conclusions? 

(1)  Systems  which  include  inlets  having  a rising  slope  of  the  curve 
total  pressure  recovery  with  increasing  volume  flow  through  the 
inlet  may  experience  both  static  and  dynamic  instability. . High- 
pressure-recovery  inlets  ingesting  boundary-layer  air  usually  will, 
have  this  characteristic  at  low  flow  rates. 

(2)  Static  instability  is  usually  observed  as  in  the  fora  of  unequal 
flow  rates  in  ostensibly  identical  branches  of  the  system.  Dynamic 
instability  occurs  in  the  form  of  regular  periodic  oscillation  of  the 
system. 

(3)  Both  static  and  dynamic  instability  are  usually  eliminated  once  all 
inlets  are  operating' at  flow  rates  well  in  excess  of  that  for  maximum 
pressure  recovery.  It  appears , however , that  instability  may  persist 
somewhat  past  the  point  of  peak  recovery  pressure  due  to  nonlinearity 
of  the  system  and  the  presence  of  finite  disturbance, 

(U)  Both  static  and  dynamic  instability  can  usually  be  eliminated  by  the 
introduction  of  throttling  devices  producing  a pressure  drop  of  the 
order  of  the  dynamic  head  of  the  entering  air. 

(5)  Although  the  static  and  dynamic  instability  have  a common  basic  cause, 
namely  the  inlet  characteristic , the  necessary  and  sufficient  condi- 
tions for  their  occurrence  are  quite  different,  and  depending  on  the 
dynamic  characteristics  of  the  system,  either,  both,  or  neither  may 
be  encountered  when  the  inlet  characteristic  is  unstable. 

(6)  Static  instability  of  the  flow  into  wide-span  slots  occurred  for 
typical  inlets  tested  in  this  program.  At  very  low  flow  rates,  the 
flow  was  reversed  at  one  end  of  the  slot  and  entering  at  much  higher 
than  average  speed.  At  the  same  time,  the  flow  was  oscillating. 

(73  The  introduction  of  splitter  plates  in  the  wide-span  slot  and  dif- 
fuser eliminated  the  dynamic  instability  but  intensified  the  static 
instability  at  low  flow  rates,  in  that  the  flow  tended  to  enter  only 
some  of  the  sections  into  which  the  duct  was  split.  When  the  total 
flow  entering  the  slot  approached  to  total  boundary-layer  quantity 
ahead  of  the  slot,  the  flow  in  the  various  sections  was  uniform. 

(8)  When  regular  oscillation  of  the  duct  system  occurred,  it  was  at  a 
frequency  generally  in  good  agreement  with  the  calculated  Helmholts 
resonator  frequency  of  the  system. 
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